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Low-mass dense cores represent the state of molecular gas associated with the earliest 
phases of low-mass star formation. Such cores are called "protostellar" or "starless," depending 
on whether they do or do not contain compact sources of luminosity. In this chapter, the first 
half of the review of low-mass dense cores, we describe the numerous inferences made about the 
nature of starless cores as a result of recent observations, since these reveal the initial conditions 
of star formation. We focus on the identification of isolated starless cores and their internal 
physical and chemical properties, including morphologies, densities, temperatures, kinematics, 
and molecular abundances. These objects display a wide range of properties since they are 
each at different points on evolutionary paths from ambient molecular cloud material to cold, 
contracting, and centrally concentrated configurations with significant molecular depletions 
and, in rare cases, enhancements. 



1. INTRODUCTION AND SCOPE OF REVIEW 

Over the last decade, dedicated observations have re- 
vealed much about the earliest phases of low-mass star for- 
mation, primarily through studies of low-mass dense cores 
with and without internal protostellar sources, e.g., proto- 
stellar and starless cores. Such cores are dense zones of 
molecular gas of relatively high density, and represent the 
physical and chemical conditions of interstellar gas just af- 
ter or prior to localized gravitational collapse. These ob- 
jects are essentially the seeds from which young stars may 
spring, and define the fundamental starting point of stellar 
evolution. Since PPIV (see Langer et al. 2000; Andre et 
al. 2000), enormous strides have been made in the observa- 
tional characterization of such cores, thanks to the increased 
interest stemming from the promising early results reported 



at PPIV and to new instrumentation. 

This review summarizes the advances made since PPIV, 
and, given the large number, it has been divided into two 
parts. In Part I (i.e., this chapter), we summarize major re- 
cent observational studies of the initial conditions of star 
formation, i.e., the individual physical and chemical prop- 
erties of starless cores, with emphasis on cores not found 
in crowded regions. This subject includes their identifi- 
cation, morphologies, densities, temperatures, molecular 
abundances, and kinematics. In recent years, the internal 
density and temperature structures of starless cores have 
been resolved, regions of chemical depletion (or enhance- 
ment) have been studied, and inward motions have been 
measured. Part II of this review (see next chapter by Ward- 
Thompson et al.) summarizes what observed individual and 
group characteristics, including magnetic field properties, 
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mass distributions and apparent lifetimes, have taught us 
about the evolution of low-mass dense cores through pro- 
tostellar formation. Despite significant progress, the obser- 
vational picture of the earliest stages of star formation re- 
mains incomplete, and new observational and experimental 
data and theoretical work will be necessary to make further 
advances. 

2. STARLESS CORES 

2.1 Definition 

Stars form within molecular gas behind large amounts of 
extinction from dust. This extinction has been used to locate 
molecular clouds (Lynds, 1962), isolated smaller clouds 
{Clemens andBarvainis, 1988), and particularly opaque re- 
gions within larger clouds (Myers et al, 1983; Lee and My- 
ers 1999). Locations within these clouds of relatively high 
density or column density, though identifiable from opti- 
cal or infrared absorption (e.g., see chapter by Alves et al), 
can also be detected using submillimeter, milhmeter, or ra- 
dio emission. For instance, spectroscopic studies of clouds 
found via extinction studies further selected those with ev- 
idence for dense gas (e.g., see Myers et al, 1983; My- 
ers and Benson 1983). These became known collectively 
as "dense cores," based primarily on whether or not they 
showed emission from NH3, indicative of densities above 
about 10* cm~^ {Benson and Myers, 1989; see also Jijina et 
al., 1999). Using IRAS data, Beichman et al. (1986) found 
that roughly half the known dense cores in clouds contained 
IRAS sources, while Yun and Clemens (1990) found that 
about 23% of the isolated smaller clouds contained infrared 
sources within; the remainder became known as "starless 
cores." Follow-up studies of clouds with sensitive submil- 
limeter continuum detection systems allowed further dis- 
crimination within the class of starless cores. Those with 
detected emission, indicating relatively high densities of 
10^~^ cm~^ (Ward-Thompson et al., 1994), were called 
"pre-protostellar" cores or later "prestellar cores." Molec- 
ular fine studies showed that prestellar cores were indeed 
more likely to show evidence for inflowing material than 
were the merely starless cores (e.g., see Gregersen and 
Evans, 2000). 

In this review, we describe the recent observational char- 
acterization of "starless cores," low-mass dense cores with- 
out compact luminous sources of any mass. (We call such 
sources "young stellar objects," regardless of whether or not 
they are stellar or substellar in mass.) By "low-mass," we 
mean cores with masses M < 10 M©, such as those found 
in the nearby Taurus or Ophiuchus clouds. Starless cores 
are important because they represent best the physical con- 
ditions of dense gas prior to star formation. Conceptually, 
we distinguish starless cores that are gravitationally bound 
as "prestellar cores" though it is very difficult to determine 
observationally at present whether a given core is bound 
or not (see below). In addition, we consider such cores 
as either "isolated" or "clustered" depending on the local 



surface density of other nearby cores and young stellar ob- 
jects. Observations of isolated cores (e.g., Clemens and 
Barvainis, 1988) are easier to interpret since they generally 
inhabit simpler, less confused envirormients, i.e., without 
nearby cores or young stellar objects. Observations of clus- 
tered cores (e.g., Motte et al., 1998) are important also be- 
cause they are more representative of star formation within 
embedded clusters, i.e., where most stars form within the 
Galaxy {Lada andLada, 2003). Finally, we describe cores 
as either "shielded" from or "exposed" to the interstellar ra- 
diation field (ISRF), depending on whether they do [e.g., 
many of those described by Myers et al, 1983] or do not 
[e.g., those cataloged by Clemens andBarvainis, 1988] re- 
side within a larger molecular cloud. 

Observational hmitations sttongly influence how dense 
cores are defined. For example, low resolution data may 
reveal cores over a wide range of mass, but higher reso- 
lution data may reveal substructure within each, with their 
own range of (lower) mass, i.e., cores within cores. In ad- 
dition, different tracers may probe different kinds of struc- 
ture within clouds, due to varying density, temperature, or 
chemistry. Also, whether or not a given dense core con- 
tains a compact infrared-luminous source is ultimately a 
matter of observational sensitivity (e.g., L1014; see below). 
Finally, learning whether or not a given core is prestellar 
(i.e., gravitationally bound) requires accurate determination 
of its physical state (mass, temperature, and internal mag- 
netic field). Masses alone can be uncertain by factors of 
^^3 due to uncertainties in distances, dust opacities, molecu- 
lar abundances, or cahbration. Various new instrumentation 
of higher resolution and sensitivity will soon be available, 
making such observational characterization more tenable. 

2.2 Identification 

Observations of starless cores have focused mostly on 
examples within 1 kpc, where the highest mass sensitivi- 
ties and linear resolutions are possible. Many starless cores 
have been detected serendipitously, adjacent to previously- 
known protostellar cores or young stellar objects (e.g., see 
Wilson etal, 1999; Williams and Myers, 1999; Swift et al, 
2005). Only recently, more systematic mapping has begun 
to identify larger numbers of starless cores in less biased 
surveys. 

Recent large line maps across significant portions of 
nearby molecular clouds using transitions that trace rela- 
tively high densities have been useful in identifying dense 
cores. For example, 179 cores in clouds within 200 pc were 
identified from substructure detected in C^*0 (1-0) maps of 
L1333 in Cassiopeia {Obayashi et al, 1998), Taurus {On- 
ishi et al, 1998), the Southern Coalsack region (Kato et 
al, 1999), tiie Pipe Nebula region {Onishi et al, 1999), 
Chamaeleon {Mizuno et al, 1999), Lupus {Hara et al, 
1999), Corona AustraHs (Yonekura et al, 1999), and Ophi- 
uchus (Tachihara et al, 2000). Tachihara et al (2002) iden- 
tified 76% of these cores as "starless" (i.e., without associ- 
ated IRAS sources; see below) with masses of 1-100 M0 
(on average, 15 M©). Using a still higher density molecu- 
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Fig. 1. — The prestellar core L1551-MC shown in NH3 (1,1) integrated intensity (greyscale) and CCS(32-2i) integrated intensity 
(contours). Positions of nearby objects are shown as symbols (see Swift et al, 2005). 



lar tracer, Hi3cO+ (1-0), Onishi et al. (2002) identified 55 
"condensations" in Taurus, 80% of which were similarly 
deemed "starless" with estimated masses of 0.4-20 Mq. 
Using line maps to identify cores can be problematic, how- 
ever, because of hne optical depths and chemical evolution 
within the cloud. For example, lines Uke C^^O (1-0) can be 
too opaque to sample most of the material in dense cores. 
Furthermore, this line and H^'^CO+ (1-0) can be affected 
by the significant depletion of the source molecules, again 
making it difficult to sample their column densities. In con- 
trast, N-bearing molecules, in particular NH3 and N2H+, 
appear to survive at higher densities than C-bearing species, 
and in fact they have been used to probe core nuclei (see 
Casein et al., 2002c; also Section 4). 

Submillimeter or millimeter continuum emission maps 
of molecular clouds can trace the high column densities 
associated with starless cores without the potentially con- 
fusing effects of Une opacity or chemical evolution. These 
studies trace the dust component, which emits radiation at 
long wavelengths even at the very low dust temperatures 
(Td) in dense cores (~10 K), though the emission is sensi- 
tive to Td. The optical properties of grains are also likely to 
evolve at relatively high densities due to ice mantle growth 
and coagulation, but these effects are much less severe than 
those for gas phase tracers (see Section 4). Recent single- 
dish telescope studies include the inner regions of Ophi- 
uchus (Motte et al, 1998; Johnstone et al, 2000; John- 
stone et al., 2004; Stanke et al, 2006), northern Orion A 
(Chini et al, 1997; Lis et al, 1998; Johnstone and Bally, 
1999), NGC 2068/2071 in Orion B {Mitchell et al, 2001; 
Johnstone et al, 2001; Motte et al, 2001), NGC 1333 in 
Perseus (Sandell and Knee, 2001), and R Corona Australis 



(Chini et al, 2003; Nutter et al, 2005). For example, Motte 
et al. (1998) found 58 "starless clumps" in their map of 
Ophiuchus with masses of 0.05-3 M©. Even larger regions 
(e.g., 8-10 deg^) have now been mapped (Enoch et al, 
2006; Young et al, 2006; Hatchell et al, 2005; Kirk et al, 
2006), providing expansive (though shallow) coverage of 
the Perseus and Ophiuchus clouds. Collectively these stud- 
ies have revealed hundreds of compact continuum objects. 
In those cases where comparisons with near- to far-infrared 
data were made, more than ^60% were deemed "starless." 
The fraction of these that are prestellar is not known, how- 
ever, and determining this (through appropriate molecular 
tracers; see Section 4) is a key goal of future research. 

In both line and continuum maps, the methods by which 
objects are identified have varied widely. Some studies have 
identified objects by eye, while others have used more ob- 
jective automated structure identification algorithms (e.g., 
Clumpfind by Williams et al, 1994), multiscale wavelet 
analyses, and even hybrid techniques. In addition, objects 
identified in continuum studies are significantly more com- 
pact in scale than those found in the line maps described 
above. First, angular resolutions of the data have substan- 
tially differed, with hne and continuum maps typically hav- 
ing > 1' and < 1' resolutions respectively. Second, the spa- 
tial sampling of the data also have substantially differed, 
though single-dish telescopes were used for both. Line 
data are obtained through frequency or position switching 
to line-free locations, allowing large-scale information to 
be retained, but continuum data are obtained through chop- 
ping (or sky background removal for instruments that do not 
chop), filtering out large-scale information (e.g., over 2-10' 
depending on the instrument.) Indeed, many continuum ob- 
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jects can be found embedded within objects identified via 
line surveys. Continuum cores have been identified, how- 
ever, using relatively low-resolution data, e.g., from IRAS 
(Jessop and Ward-Thompson, 2000) or ISO {Tdth et al, 
2004). 

Identification of objects found through fine or contin- 
uum observations as starless cores is ultimately dependent 
upon the availability and sensitivity of ancillary data at the 
time of analysis. For example. Figure 1 shows the starless 
core L1551-MC in Taurus recently discovered by Swift et 
al. (2005) in NH3 and CCS emission, with the positions 
of nearby, non-coincident objects hsted in the SIMBAD 
database (as of late 2004). Far-infrared data from IRAS 
were used traditionally to find evidence for embedded pro- 
tostars associated with detected cores, but more recently 
wide-field near-infrared data, including those from 2MASS, 
have been also used. Such wide-field data can be insuf- 
ficiently sensitive to detect the youngest, most-embedded 
protostars, however, given high extinctions (e.g.. Ay > 
30). Near-infrared observations targeting specific cores can 
attain much higher sensitivities but these have not been 
done extensively. Some targeted studies, e.g., Allen et al. 
(2002) or Murphy and Myers (2003), however, did not de- 
tect low luminosity protostars in various cores despite rel- 
atively high sensitivity. The extraordinarily high sensitiv- 
ity of the Spitzer Space Telescope in the mid-infrared has 
revealed about 12 low luminosity protostars within cores 
previously identified as starless, e.g., in L1014 by Young et 
al. (2004a) and in LI 148 by Kaujfinann et al. (2005; (see 
also Bourke et al, in preparation; Huard et al, in prepara- 
tion). These are candidates to join the class of Very Low 
Luminosity Objects (VeLLOs), sources with L,,i/ < 0.1 L© 
that are located within dense cores. (Note that V^mt is the 
luminosity of the object, in excess of that supplied by the 
ISRF.) Based on this definition, the previously identified, 
extremely young, low-luminosity Class protostar IRAM 
04191 {Andri etal, 1999) is also a VeLLO {Dunham etal, 
in preparation). 

In addition to near- to far-infrared imaging, deeply em- 
bedded protostars may be identified using other signposts 
of star formation not affected by high extinction, includ- 
ing compact molecular outflows, "hot core" molecular line 
emission, masers, or compact thermal radio emission. For 
example, Yun et al. (1996) and Harvey et al. (2002) used the 
latter method to hmit luminosities of protostellar sources 
within several starless cores using limits on radio emission 
expected from shock-ionized gas, i.e., from outflows. 

3. DENSITY AND TEMPERATURE STRUCTURES 

3.1 Probing Physical Structure 

In this section, we describe the morphologies and in- 
ternal configurations of density and temperature of starless 
cores as derived in recent studies. Much progress in un- 
derstanding dense core structure has come from analyses of 
far-infrared to millimeter continuum emission observations. 



typically of relatively isolated and exposed examples. Con- 
tinuum data have been used preferentially over molecular 
fine data because interpretation of the latter can be com- 
plicated by large internal molecular abundance variations 
(see Section 4). Near- or mid-infrared absorption of back- 
ground emission by cores can be also a powerful probe of 
the density structure since extinction is independent of tem- 
perature and can trace those locations where dust colunm 
densities are too low to be easily detected in emission at 
present (e.g., see Bacmann etal., 2000). Results about core 
structure derived from "extinction mapping" have agreed 
generally with those from emission studies, with some ex- 
ceptions (see below and the chapter by Alves et al.) 

The observed emission from cores depends on the entan- 
gled effects of temperature, density, properties of the tracer 
being used, and the observational technique. For molecu- 
lar lines, the abundance and excitation of the tracer varies 
within the core. For continuum emission from dust, the in- 
terpretation depends on dust properties and the characteris- 
tics of the ISRF. Two general approaches are taken: in the 
first, simphfying assumptions are made and the observa- 
tions are inverted to derive core properties, usually averages 
over the line of sight and the beam; in the second, specific 
physical models are assumed and a self-consistent calcu- 
lation of predicted observations is compared to the actual 
observations, allowing the best model to be chosen. Both 
approaches have their strengths and weaknesses. 

Self-consistent models must compute the temperature 
distribution for a given density distribution, dust proper- 
ties, and radiation field. For starless cores, only the ISRF 
heats the dust. In their models, Evans et al. (2001) included 
the "Black-Draine" (BD) ISRF, one based on the COBE re- 
sults of Black{l99A), which includes the cosmic microwave 
background (CMB) blackbody and the ISRF in the ultravi- 
olet of Draine (1978; cf. van Dishoeck 1988), rather than 
flie older, pre-COBE model of Mathis et al. (1983; MMP). 
The BD and MMP models can differ significantly, e.g., up 
to a factor of 13 between 5 /im and 400 /im. For dust opac- 
ities, those currently in use have been theoretically derived 
from optical constants of grain and ice materials and in- 
clude features such as mid-infrared vibrational bands (e.g., 
9.7 iim Si-O stretch) and tend toward a power-law decrease 
of opacity at submilUmeter wavelengths (k^ oc A~'^). Fig- 
ure 2a shows five dust opacity models: OH5 and OH8 are 
models of grains that have coagulated for 10^ yr and ac- 
quired varied depths of ice mantles (Ossenkopf and Hen- 
ning, 1994); "Draine-Lee" are the "standard" interstellar 
medium dust opacities derived from silicate and graphite 
grains by Draine and Lee (1984); MMP opacities are de- 
rived from an empirical fit to observed dust properties in 
star-forming regions; and Pollack opacities are derived from 
an alternative grain composition based on siUcate grains 
(olivine and orthopyroxene), iron compounds (troilite and 
metallic iron), and various organic C compounds {Pollack 
et al., 1994). The mass opacity {Kv) at long wavelengths 
can vary by one order of magnitude between the opacity 
models. 



4 





»t (AU) 



Fig. 2. — a) Opacity vs. wavelength from selected dust models (see Section 3.3; Shirley et al, 2005). Values (ki^) and spectral indices 
of opacities at 850 fim are tabulated. Units of k^, are cm^ and the p at 850 fim were obtained from linear regressions of 
respective model values between 350 /um and 1.3 mm. b) The ratio of average density (n) to density at the BBS boundary radius (rib) vs. 
BBS boundary radius (Rb) for families of BBS of given ric {Shirley and Jorgensen, in preparation). The solid horizontal line corresponds 
to the stability criterion of Lombardi and Berlin (2001). 



These ingredients then provide inputs to radiative trans- 
port codes (e.g., Egan et al, 1988; Ivezic et al.,, 1999), 
which attempt to fit simultaneously the observed intensity 
profiles, Iu{t) (most sensitive to variation in n(r)), and 
spectral energy distribution (SED, most sensitive to the total 
mass and the assumed dust opacity). To test various ISRF 
and dust opacity models for the L1498 starless core, Shirley 
et al. (2005) varied the BD ISRF by a scaling factor, .s,.,, /, 
and tested eight dust opacity models, including those shown 
in Figure 2a. In their models, the resultant SED was more 
sensitive to the ISRF than the resultant normalized contin- 
uum intensity profiles, and they found best fits with Sisrf 
= 0.5-1.0. Furthermore, only the OH5 and OH8 opacity 
models were consistent with the observed SED. 

Chemical differentiation observed in starless cores (see 
Section 4) should lead to variations in ice compositions, and 
hence dust opacities. Indeed, Kramer et al. (2003) found 
in IC 5 146 an increase in dust emissivity by a factor of 4 
between 20 K and 12 K. To date, variations of dust opac- 
ity with radius have not been self-consistently included in 
radiative transfer models. Moreover, most radiative trans- 
fer models are ID and are constrained by ID intensity pro- 
files, although starless cores can be quite ellipsoidal in 3D 
shape (see §3.2). If magnetic fields are dynamically impor- 
tant, they require core geometries that are inherently multi- 
dimensional since the magnetic field only provides support 
perpendicular to the direction of the field fines (e.g., Li and 
Shu, 1996). Promising work on modeUng cores with 3D 
radiative transfer codes has begun (e.g.. Doty et al., 2005a, 
2005b; Steinacker et al., 2005). A reanalysis of L1544 core 
data by Doty et al. (2005a) using a 3D code yielded results 



that were overall similar to those found from ID models but 
with reduced size and mass. 

3.2 Morphology and Density Structure 

Starless cores vary widely in morphology or shape, rang- 
ing between filamentary to compact and round. Cores are 
seen as 2D projections against the sky, making it difficult to 
constrain their 3D shape. Furthermore, core morphologies 
can depend on the wavelength and the angular resolution of 
observation (e.g., see Stamatellos et al, 2004). Observed 
morphologies also depend upon the intensity level chosen 
as the boundary between cores and their surroundings, and 
this can be affected strongly by observational sensitivity. In 
general, cores that are round (in projection) have been typ- 
ically interpreted as being spherical in shape, while cores 
that are elongated (in projection) have been typically inter- 
preted as being either oblate or prolate spheroids. Statistical 
analyses of core map aspect ratios indicate that most iso- 
lated cores are prolate rather than oblate {Myers etal, 1991; 
Ryden, 1996), and many cores in Taurus share the projected 
ahgnment of the filamentary clouds in which they reside 
{Hartmann, 2002). Jones et al. (2001), using a sample of 
cores identified from the extensive NH3 catalog of Jijina et 
al. (1999), posited that elongated cores actually have intrin- 
sically triaxial shapes that are neither purely oblate nor pro- 
late. Goodwin et al. (2002), using the same catalog, found 
starless cores have more extreme axial ratios than protostel- 
lar cores. 

Column densities within starless cores can be traced us- 
ing the optically thin millimeter or submillimeter contin- 
uum specific intensity. Derivation of absolute column den- 
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sities and conversion to volume densities requires knowl- 
edge or assumptions about dust temperatures, opacities, 
gas-to-dust ratios, geometry, and telescope beam pattern. 
Azimuthally averaged (or elliptically averaged in cases of 
extreme axial ratio), ID intensity profiles of starless cores 
are found to have "flat" slopes at radii < 3000-7500 AU 
and "steep" slopes at larger radii (see Ward-Thompson et 
at, 1994; Andre et al., 1996; Ward-Thompson et al., 1999; 
Shirley et al, 2000; Kirk et al., 2005). Such profiles depart 
from the n(r) oc distribution of the "singular isother- 
mal sphere" postulated by Shu (1977) as the initial state of 
isolated dense cores prior to gravitational collapse. Sim- 
ilar profiles have been seen also in studies of dense core 
structure utilizing the absorption of background emission 
by core dust (e.g., Alves et al., 2001; see chapter by Alves 
et al.) The observed configurations are better described by 
"Bonnor-Ebert spheres" (BES; Emden, 1907; Ebert, 1955; 
Bonnor, 1956; Chandrasekhar, 1957), non-singular solu- 
tions to the equations of hydrostatic equilibrium that can be 
critically stable in the presence of external pressure. In the 
ID isothermal approximation, the family of BES solutions 
are parameterized by the central density, Uc, and are char- 
acterized by density profiles with two regimes: a "plateau" 
of slowly decreasing density at small radii and a power-law 
decrease (~ r~^) at large radii {Chandrasekhar and Wares, 
1949). The size of the plateau scales inversely with the cen- 
tral density, i.e., R{\nc) oc \/T/nc. For reference. Figure 
2b shows the average density contrast, where n is the 
average density and Ub is the density at the BES (outer) 
boundary radius, versus BES boundary radius, Rb, for fam- 
ilies of BES of given ric {Shirley and Jorgensen, in prepara- 
tion). Thermally supported density configurations, regard- 
less of geometry, are unstable if they have ^ > 2.46457 
{Lombardi and Bertin, 2001). Such configurations may 
remain stable if non-thermal internal pressures (e.g., from 
magnetic fields or turbulence) are relevant, however; see 
Section 5, as well as the next chapter by Ward-Thompson et 
al. 

Detailed modeling of continuum intensity profiles have 
confirmed that BES are a good approximation to the den- 
sity structures of starless cores (e.g., Ward-Thompson et 
al, 1999; Evans et al., 2001; Kirk et al, 2005; Schnee 
and Goodman, 2005, but see Harvey et al., 2003). Re- 
cently, Shirley and .Jorgensen (in preparation) modeled sub- 
millimeter continuum emission from 33 isolated starless 
cores, the largest sample so far, and found a log-normal 
distribution of central densities with median logn,; = 5.3 
cm~^. This large median central density indicates that non- 
thermal support (e.g., magnetic pressure) has stabilized the 
cores, that the density scale is inaccurate due to the dust 
opacity assumed (see below), or that these cores are not 
in equilibrium. Unfortunately, it is not currently possible 
to distinguish between collapsing and static BES with dust 
emission or absorption profiles because the density struc- 
ture, for a given ric, does not significantly vary until late 
in the collapse history {Myers, 2005). Molecular probes of 



the kinematics are much better suited to resolving this issue 
(e.g., see Section 5). Nevertheless, increasing central con- 
centration within starless cores could indicate an evolution- 
ary path. Note, however, that while BES are consistent with 
current observations, this does not prove that given physi- 
cal structures are BES (cf. Tafalla et al., 2002; Ballesteros- 
Paredes etal, 2003, but see Section 5.2). 

3.3 Temperature Structures 

Early calculations of dust temperature (To) in cores 
heated only by the ISRF indicated To ~ 10 K {Leung, 
1975), cooler on the inside and warmer on the outside. Such 
a gradient in Td is unavoidable theoretically because the 
cores are optically thin to their own radiation (Apeak ~ 200 
/um) and are primarily heated externally by a local ISRF 
for which short wavelength radiation is significantly obnu- 
bilated in the center of the core. Alternative heating mech- 
anisms for dust, including cosmic rays (i.e., direct heating 
of grains, secondary UV heating from excitation of H2, and 
heating of gas followed by gas-grain heating) are not signif- 
icant compared to heating from the ISRF. More recent cal- 
culations using current ideas for the ISRF and dust opacites 
confirm the earlier work (e.g., Zucconi et al., 2001; Evans 
et al., 2001), finding To falling from about 12 K on the 
exposed surface to as low as 7 K at very small radii. Such 
low temperatures would suppress radiation even at millime- 
ter wavelengths from the center, possibly affecting interpre- 
tation of the central densities. Gradients of found in 
starless cores, however, do not have a strong effect on over- 
all core equilibria, and density profiles very different from 
isothermal BES are not expected {Evans et al, 2001; Galli 
etal, 2002). 

The temperature gradient within a core depends on the 
strength and spectral properties of the ISRF, which are af- 
fected by surrounding extinction, and the density structure 
in the core. Shielded cores have lower temperatures and 
shallower gradients than exposed cores, because the ISRF is 
attenuated and reddened {Stamatellos and Whitworth, 2003; 
Andre et al, 2003; Young et al, 2004b). The reddening 
leaves more of the heating to the longer wavelength pho- 
tons, which "cook" more evenly. Cores of relatively low 
central concentration or low central density should have 
shallower temperature gradients due to lower extinction 
{Shirley et al, 2005). Calculations in higher dimensions 
also show smaller gradients. For clumpy cores, the heat- 
ing is dominated by unobscured lines of sight {Doty et al, 
2005b). Such clumpy cores can have cold patches {To < 6 
K) within an overall flatter distribution of Tq. 

Observational tests do indicate warmer Tp on the out- 
side of cores. Langer and Willacy (2001), Ward-Thompson 
et al (2002) and Pagani et al. (2004, see also 2003), each 
using ISO data, found evidence for To gradients in exposed 
cores, with innermost To < 10 K. Limited spatial resolu- 
tion at far-infrared wavelengths, however, does not allow 
direct tests of the predictions of very low To at very small 
radii. One indirect test comes from studies of starless cores 
in extinction (e.g., Bacmann et al, 2000), which are tem- 
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perature independent. These studies also show "flat" inner 
density gradients, and hence these gradients are not due to 
decreased To in the centers. Another possible indirect test 
is the high-resolution study of molecular lines (see below). 

Examining gas kinetic temperature (Tk), Goldreich and 
Kwan (1974) showed it is forced to be very close to Td 
at high densities (n > 10® cm~^). This coupling fails for 
less dense cores and for the outer parts of even denser cores 
{Goldsmith, 2001; Galli et al, 2002). At first, Tk may 
drop below To, but cosmic ray heating prevents Tk from 
dropping too low. If the core is not very deeply shielded, 
however, photoelectric heating causes Tk > To in the outer 
layers (e.g.. Young et al., 2004b). 

Reductions in the ISRF affect the chemistry in the outer 
layers {Lee et al, 2004b) and change the gas kinetic tem- 
perature {Tk) profile as well by decreasing photoelectric 
heating. Conversely, cores in regions of high ISRF will 
have warmer surfaces in both dust and gas {Lis et al, 2001; 
J0rgensen et al., 2006). Since this heating is caused by 
ultraviolet light, it is very sensitive to small amounts of 
shielding. CO is a particularly good constraint on the short- 
wavelength end of the ISRF (e.g., Evans et al, 2005), but, 
with chemical models (Section 4) other species can also 
be excellent probes (e.g.. Young et al, 2004b; Lee et al, 
2004b). 

Observationally, Tk can be determined directly from 
CO transitions, which probe the outer layers, and by in- 
version transitions of NH3, a good tracer of dense core ma- 
terial (e.g., see Jijina et al, 1999 and references therein), 
which probe deeper layers. Analyses of such lines from 
starless cores have found nearly constant values of Tk » 
10 K (e.g., Hotzel et al, 2002 for B68 and Tafalla et al, 
2004 for L1498 and L1517B). If Tk = Td, and Td de- 
creases toward the center, one might expect to see changes 
in Tk- For L1498 and L15 17B, however, Galli et al (2002) 
found Tk was ~2 K higher than Td for cores with cen- 
tral densities ~10® cm'^ (i.e., Uke L1498 and L1517B), 
and that neither Tk nor Td varied significantly within the 
inner core radii. Shirley et al. (2005), however, found a 
much lower central density (i.e., nc = 1 — 3 x 10* cm~^) 
for L1498, which would explain the difference of Tk and 
Td. Also, recent VLA observations of NH3 (1,1) and (2,2) 
across the L1544 prestellar core by Crapsi et al (in prepa- 
ration) have revealed a clear Tk gradient that is consistent 
with the dust temperature profile derived from models of 
submilUmeter continuum emission with a central density of 
3 X lO'^ cm-3. 

4. CHEMICAL CHARACTERISTICS 

Structures of starless cores in principle can be probed 
using molecular emission lines, typically rotational tran- 
sitions excited at the relatively low temperatures and high 
densities of such objects. Indeed, line data can comple- 
ment continuum data to provide new constraints to models 
(e.g., Jessop and Ward-Thompson, 2001). Line studies also 



have the advantage of tracing kinematic behavior of starless 
cores. In practice, however, line observations are affected 
by the abundance variations within the cores and by the op- 
tical depths of the specific lines used; unlike submillime- 
ter or millimeter continuum emission, the optical depths of 
molecular lines are often quite large. To circumvent this 
problem, Unes from rare isotopologues are observed and 
interpreted using non-LTE radiative transfer codes, which 
take into account the core chemical structure. 

Although the basic chemical processes in clouds have 
been known for some time (e.g., formation of H2 on dust 
grains {Gould and Salpeter, 1963; Jura, 1975), ioniza- 
tion of H2 by cosmic rays {Solomon and Werner, 1971), 
ion-molecule reactions {Herbst and Klemperer, 1973), and 
grain-surface chemistry {Watson and Salpeter, 1972; Allen 
and Robinson, 1977)), there are still several uncertain- 
ties for denser regions, one of which is grain-surface pro- 
cesses. Theoretical studies predicted molecular depletion 
onto grain surfaces (e.g., Leger, 1983). While adsorbed 
molecules have long been observed in infrared absorption 
bands (e.g., Gillett and Forrest, 1973; see also Pontoppi- 
dan et al, 2005), only in the past few years have molecu- 
lar line observations found clear and quantitative evidence 
for gas-phase depletion (e.g., Kuiper et al, 1996; Willacy 
et al, 1998; Kramer et al, 1998, see Section 4.1). En- 
hanced deuterium fractionation, which should accompany 
molecular depletion (e.g.. Brown and Millar, 1989), is also 
found in both starless and protostellar cores (see Section 
4.2). Uncertainties still remain, however, and further up- 
dates of chemical models and detailed comparisons with 
new observational data are needed. 

4.1 Molecular Freeze-out 

The dominant gas phase constituent in molecular clouds 
is H2 but its line emission is not used to trace the interi- 
ors of molecular clouds and cores. H2, as a homonuclear 
molecule, has no electric dipole moment and thus has no 
dipole transitions between rotational states. Quadrupole 
emission between states is possible but the upper level of 
the first allowed transition is 5 12 K above ground (since H2 
is a light molecule) and is only excited where gas is suitably 
hot, e.g., on cloud surfaces and not in cloud interiors. 

CO and its isotopologues ^^CO, C^^O or C^'^O are com- 
monly used as surrogate tracers to H2 since they are rel- 
atively abundant and intermixed with H2, and have rota- 
tional transitions excited at the ambient densities in molec- 
ular clouds, i.e., lO^^"* cm^"*. Recent studies, however, 
have found significant depletions of CO abundance within 
the innermost regions of starless cores by comparing Une 
observations of CO isotopologues with the submillimeter 
dust continuum or extinction maps of cores. For example, 
Bacmann et al (2002) found evidence for CO depletions by 
factors of 4-15 in a sample of 7 cores. Furthermore, Caselli 
et al (1999) and Bergin et al (2002), using C^^O and C^^O 
respectively, found CO abundance depletions by factors of 
--1000 and 100 in the centers of the LI 544 and B68 cores 
respectively, i.e., where n > 10® cm~^. 
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The most common explanation for these abundance vari- 
ations is that CO and its isotopologues adsorb, i.e., "freeze- 
out," easily onto grain mantles at high densities and dust 
temperatures < 20 K. Other molecules are also similarly 
depleted onto dust grains within low-mass cores. For ex- 
ample, Bergin et al. (2001), Young et al. (2004b), Lai et al. 
(2003) and Ohashi et al. (1999) found observational evi- 
dence for depletions of CS, H2CO, C3H2 and CCS respec- 
tively toward various dense cores. More recently, Tafalla et 
al. (in preparation) carried out a multimolecular study of the 
L1498 and L15 17B cores and found evidence for the deple- 
tion of 1 1 species, including all the aforementioned species 
plus DCO+, HCN, HC3N, HCO+, CH3OH and SO. 

Water is also clearly frozen onto dust grains. First, solid 
water is seen (with abundances ^10^^ relative to H2) in 
absorption within the near-infrared spectra of stars located 
behind molecular clouds and those of embedded protostars 
(e.g., van Dishoeck, 2004 and references therein). Indeed, 
water ice is the major compound of grain mantles, followed 
by CO (e.g., Chiar et al., 1995; Ehrenfreimd and Cham- 
ley, 2000). Second, the recent SWAS mission furnished 
stringent upper Umits to abundances of water vapor (< 1 0~^ 
relative to H2; Bergin and Snell 2002) that are orders of 
magnitude lower than values predicted by gas-phase chem- 
ical models that did not account for the gas-grain interac- 
tion and the consequent molecular freeze-out. Other ices on 
grain surfaces recently detected through absorption spectra 
of background stars include CO2 (Bergin et al., 2005), and 
probably HCOOH and NH+ (Knez et al., 2005). 

The degree to which a given molecule depletes onto 
grains is Ukely a function of its respective binding en- 
ergy onto grain surfaces, but other chemical factors may be 
also important. For example, carbon-chain molecules such 
as CCS and C3H2, which are called "early-time" species 
because in gas-phase chemical models they form before 
atomic carbon is mainly locked in CO (Herbst and Leung, 
1989), can decrease in abundance due to gas-phase reac- 
tions, especially with atomic oxygen, which tend to convert 
all C-bearing species into CO. When CO and the reactive O 
start to freeze-out, carbon-chain species increase their abun- 
dance again, showing a "late-time" secondary peak, limited 
in time by their own adsorption onto grains (Ruffle et al, 
1997, 1999; Li etal, 2002; Lee etal, 2003). 

Nitrogen-bearing molecules, e.g, NH3 and N2H+, have 
been shown to be quite resihent tracers of dense cores, rel- 
ative to carbon-bearing molecules. They are called "late- 
time" species; they both have N2 as a parent molecule, 
which forms slowly (relative to CO) in the gas phase via 
neutral-neutral reactions. As described in Section 2, NH3 
transitions have long been tracers of dense cores since these 
can be excited at densities w 10^ cm^-^. In recent years, 
however, N2H+ transitions have grown in popularity for 
tracing dense cores with surveys by Benson et al. (1998), 
Casein et al. (2002c) and Tatematsu et al. (2004). The in- 
tensity distribution of N2H+ 1-0 in a given core closely 
matches that of the millimeter or submillimeter continuum 
(e.g.. Casein et al, 2002a) and its hyperfine structure can 



be fit to determine excitation temperatures and line opac- 
ities (under the assumption that all the hyperfine compo- 
nents have the same excitation temperature; see Womack et 
al. 1992; Caselli et al. 1995). Within 5 cores, Tafalla et 
al. (2002; see also Tafalla et al, 2004) found inner zones of 
CO and CS depletion by factors > 100 but NH3 abundances 
that increased by factors of ~1-10 towards the core centers 
and constant N2H+ abundances. 

The relatively high abundances of NH3 and N2H"'" in 
starless cores were thought to be due to the low binding en- 
ergy of their parent N2 relative to CO (Bergin and Langer, 
1997), but recent laboratory experiments by Oberg et al. 
(2005; see also Bisschop et al., 2006) have found the bind- 
ing energies to be quite similar Another possible explana- 
tion is a slow formation of N2 molecule in the gas phase 
(Caselli andAikawa, in preparation). In the initial forma- 
tion of molecular clouds, the N atom would be the dominant 
form of nitrogen. The nitrogen atom has a lower adsorp- 
tion energy than N2 and CO, and thus is more resilient to 
adsorption. Formation of N2 is much slower than the CO 
formation in the gas phase, and hence many N atoms are 
still available when the adsorption starts (i.e., when the core 
age exceeds the collisional timescale between gas and dust). 
Adsorption of N2 is compensated by the formation in the 
gas phase. The high abundance of their mother molecule 
helps N2H+ and NH3 to maintain their abundances. Two 
additional factors lead to high N2H+ abundance. First, 
N2H+ does not directly adsorb onto grains; it recombines 
on charged grain surfaces and returns to the gas-phase as 
N2. Second, the main reactants of N2H+, e.g., CO and 
electrons, decrease as the core density rises. Indeed, in 
Model A of Aikawa et al. (2001), depletion of N2H+ and 
NH3 is much less significant than that of CO, even though 
the adsorption energies of CO and N2 are assumed to be 
equal. The central NH3 enhancement relative to N2H+ can 
be caused by CO depletion. In the outer regions with high 
CO abundance, N2H+ mainly reacts with CO to produce 
N2, while, in the central regions with heavy CO depletion, 
it recombines to produce NH as well as N2 (Geppert et al, 
2004), which is transformed to NH3 (Aikawa etal., 2005). 

Nitrogen-bearing species are not immune to depletion 
and even N2 will be eventually adsorbed onto grains at high 
densities in starless cores (e.g., Bergin and Langer, 1997; 
Aikawa et al., 2001). Bergin et al. (2002) and Crapsi et al. 
(2004) found evidence for small zones of N2H+ depletion 
(by factors of ~2) coincident with the centers of the ex- 
posed cores in B68 and LI 52 IF, respectively. In addition, 
Di Francesco et al. (2004) and Pagani et al. (2005) found 
evidence for N2H+ depletion in the Oph A and L183 cores, 
respectively. Belloche and Andre (2004) and Jprgensen et 
al. (2004) found N2H+ depletion within cold IRAM 04191 
and NGC 1333 IRAS 2 protostellar cores respectively, but 
outflows may have also partially affected the local chem- 
istry if CO has been released from grain mantles near the 
protostars. Further, Caselli et al. (2003) and Walmsley et 
al. (2004) have suggested that all molecules with C, N or O 
can be utterly depleted in the innermost regions of denser 
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Fig. 3. — A schematic representation of eventual molecular differentiation within a starless core. The external shell of the core (where 
n(H2) ~ 10* cm~^) can be traced by CO, CS and other carbon bearing species. At radii < 7000 AU, where n(H2) ~ 10* cm~^, CO 
and CS disappears from the gas phase and the best gas tracers are NH3 and N2H"''. At higher densities, deuterated species become quite 
abundant and, when the density exceeds ^lO'' cm~'^, the chemistry will be dominated by light molecular ions, in particular and its 
deuterated forms, as well as H+ (e.g., Walmsley et al, 2004). 



cores based on interpretation of the H2D+ emission in the 
L1544 core. Indeed, molecules without heavy elements, 
e.g., H2D+, may be the only remaining molecular tracers 
of such regions. 

Figure 3 schematically summarizes how the eventual 
molecular differentiation within a starless core may look. 
At radii of -7000-15000 AU, CO and CS are still mainly 
in the gas phase and the main molecular ion is HCO+, with 
which one can deduce a stringent lower limit of the electron 
density (e.g., Caselli et al., 1998; see Section 4.3). The frac- 
tion of atomic carbon in the gas is still large, as testified by 
the large observed abundances of carbon-chain molecules, 
such as CCS (e.g., Ohashi et al, 1999). Deeper in the 
core (~5000-7000 AU), where the density approaches val- 
ues of the order of 10^ cm~^, CO and CS disappear from 
the gas-phase because of the freeze-out onto dust grains. 
The physical and chemical properties (as well as kinemat- 
ics) are better traced by N-bearing species, in particular 
NH3. Within the central 5000 AU, deuterium fractionation 
takes over (see Section 4.2), and N2D+ becomes the best 
probe {Casein et al, 2002a). NH3 is still abundant in these 
regions, however, as suggested by the observed increase 
of the NH3 abundance toward core centers {Tafalla et al, 
2002). At r < 2500 AU, where n(H2) > 10*^ cm-^, all neu- 
tral species are expected to freeze-out in short time scales 
(<1000 yr) and light species, such as H3 and its deuterated 
forms are thought to dominate the chemistry and the degree 
of ionization (Caselli et al, 2003; Vastel et al, 2004). 

The schematic picture of Figure 3 of course depends on 
the time spent by a starless core in this condensed phase, so 



that one expects to find less significant depletion and more 
typical cloud chemistry in those objects that just entered 
this phase (e.g., possibly the L1521E core; see Tafalla and 
Santiago, 2004). Moreover, recent VLA observations of 
NH3 (1,1) and (2,2) by Crapsi et al. (in preparation) have 
shown that NH3 is still present in the central 800 AU of 
the LI 544 core, where the gas density is a few times 10^ 
cm^"^. If no efficient desorption mechanisms are at work, 
the time spent by the L1544 core nucleus in its high density 
phase may be < 500 yr (see Section 4.4). More examples 
are needed to refine chemical models of starless cores. 

4.2 Deuterium Fractionation 

Another important chemical process recently identified 
within starless cores is deuterium fractionation, i.e., the 
enhancement of deuterated isotopologues beyond levels 
expected from the elemental D/H ratio of —1.5 x 10~^ 
(Oliveira et al, 2003). For example, in a sample of dense 
cores, Bacmann et al (2003) found a D2CO/H2CO column 
density ratio between 0.01 and 0. 1 . In addition, Crapsi et al. 
(2005) found N2D+/N2H+ ratios between 0.05 and 0.4 in a 
similar sample of cores. Deuterium fractionation is related 
to core temperature and CO depletion (e.g., Dalgarno and 
Lepp, 1984). Species such as and CH^ are enriched 
in deuterium in cold clouds because of the difference in 
zero-point energies between deuterated and non-deuterated 
species and rapid exchange reactions such as + HD — > 
H2D+ + H2 (e.g., Millar et al, 1989). The enrichments are 
propagated to other molecules by chemical reactions. At 
high densities, heavy element molecules like CO will also 
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deplete onto grains, reducing the destruction rate of H2D+ 
and further increasing the H2D+/Hj|^ ratio. For example, 
the N2D+/N2H+ ratio is higher than the DCO+/HCO+ ra- 
tio by a factor of about 5 in the L1544 core (Caselli et al., 
2002b) because N2H+ can trace the central region with CO 
depletion, while HCO+, which is produced by the protona- 
tion of CO, is not abundant at the core center. 

Multiply deuterated species, such as HD|^ and D;^, can 
be similarly increased in abundance (Roberts et al., 2003; 
see the chapter by Ceccarelli et al.). For example, the 
LI 544 and IRAS 16293E core centers are also rich in H2D+ 
and HD2 , as confirmed respectively by Caselli et al. (2003) 
and Vasteletal. (2004). 

The sensitivity of the D/H ratio to physical conditions 
and molecular depletion could make it a potentially useful 
probe of chemical or dynamic evolution of cores. For ex- 
ample, Crapsi et al. (2005) observed 14 cores and found 
that deuterium enhancement was higher in those that are 
more centrally concentrated and have larger peak H2 and 
N2H+ column densities, i.e., arguably those more likely to 
collapse into stars. 

4.3 Fractional Ionization 

Observation of molecular ions, combined with chemi- 
cal models, are used frequendy to derive ionization degree, 
Xe- This value is important for dynamics since the amounts 
of free charge within cores determine the relative influence 
on core evolution of ambipolar diffusion, the gravitational 
inward motion of neutral species retarded by interactions 
with ionic species bound to a strong magnetic field (e.g., 
see Ciolek and Mouschovias, 1995). For example, Caselli 
etal. (1998; see also Williams et al, 1998) found, for dense 
cores observed in DCO"*" and HCO+, values between 
10~^ and 10~^, with associated ambipolar-diffusion/free- 
fall timescale (AD/FF) ratios of 3-200 and no significant 
differences between starless and protostellar cores. In the 
L1544 core, using N2D+ and N2H+, Caselli et al. (2002b) 
found a significantly smaller x^, ^10"^, and an associated 
AD/FF ratio of ~1, suggesting this particular core is near 
dynamical collapse. 

Estimates of x^, however, depend on many assumptions, 
e.g., the cosmic-ray ionization rate and molecular depletion 
degrees, as well as the accuracy of the chemical models 
used. In particular, the cosmic-ray ionization rate, C,, typ- 
ically assumed to be around 1-3 x 10~^^ s~^ in molec- 
ular clouds (e.g., Herbst and Klemperer, 1975; van der 
Tak and van Dishoeck, 2000), seems to be closer to 10~^^ 
s~^ in diffuse clouds (McCall et al, 2003). Therefore, a 
large variation of ( is expected in the transition between dif- 
fuse and molecular clouds (e.g., Padoan and Scalo, 2005) 
and, maybe, between dense cores with different amounts 
of shielding. Such variation may, in turn, cause signifi- 
cant variations of the ionization degree between cores (e.g., 
Caselli et al, 1998) with consequent differences in AD/FF 
time scales. Padoan et al. (2004) have further argued that 
the ISRF strengths (i.e., ^v) and core age variations may 
account for at least some of the large range of x^ seen 



in dense cores. Walmsley et al. (2004) and Flower et al. 
(2005) have also suggested that the charge and size distri- 
butions of dust grains impact ionic abundances, and thus 
the local AD/FF ratio, in the centers of cores depleted of 
heavy-element species. 

4.4 "Chemodynamical" Evolution 

Depletions or enhancements of various molecular species 
in starless cores can be used as probes of dynamical evo- 
lution, since chemical and dynamical timescales can sig- 
nificantly differ. For example, if the dynamical timescale 
of a core is shorter than the adsorption timescale of CO 
(^104( i"°['g^p )(Mg)0.5 yr)^ it not experience CO 

depletion. Recent theoretical studies have investigated the 
evolution of molecular abundance distributions in cores 
(e.g., in L1544) using models including Larson-Penston 
collapse, ambipolar diffusion of magnetized cores, and con- 
traction of a critical BBS (Aikawa et al, 2001, 2003, 2005; 
Li et al, 2002; Shematovich et al, 2003; Lee et al, 2004b), 
showing that chemistry can be used to probe, although not 
uniquely, how the core contracts. 

It is important to note that not all starless cores display 
similar chemical differentiation. Hirota et al (2002) found 
that CCS is centrally peaked in the L1521E core, while it 
is depleted at the center of the L1498 core (Kuiper et al, 
1996), although the central density of the L1521E core (^3 
X 10^ cm~^) is higher than that of the L1498 core (~1 
X 10^~^ cm~^; see Shirley et al., 2005; Shirley and Jor- 
gensen, in preparation; Tafalla et al, 2002, 2004, 2005). 
The L1521E core also has low NH3 and N2H+ abundance, 
low DNC/HNC ratio, and no CO depletion (Hirota et al., 
2001, 2002; Tafalla and Santiago, 2004). The core is ap- 
parently "chemically young"; these molecular abundances 
are reproduced either at early stages of chemical evolu- 
tion in a pseudo-time dependent model (i.e., assuming that 
the density does not vary with time), or in a dynamical 
model where the contraction or formation timescale of a 
core is smaller than the chemical timescale (e.g., Lee et 
al, 2003; Aikawa et al, 2001; 2005). Other candidates 
of such "chemically young" cores are those in L1689B, 
L1495B, and L1521B (Lee et al, 2003; Hirota et al, 2004). 
Further determinations of kinematics and physical condi- 
tions in these cores and "chemically old" cores (e.g., that in 
L1498) will be very useful toward understanding chemical 
evolution in starless cores. 

5. BULK AND TURBULENT MOTIONS 

5.1 Bulk Motions 

Molecular line profiles can be effective tracers of starless 
core kinematics, provided appropriate lines are observed 
and modeled. For example, internal velocity gradients in 
dense cores can be derived from the variation of Une pro- 
files across cores. Recent interpretations of such gradients 
as rotational in origin include those made by Barranco and 
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Fig. 4. — Spectra of the F1,F = 2,3— >1,2 hj^perfine component of N2H+(l-0) (histograms) from the starless cores in L694-2 (left) 
and L1544 (right) from Williams et al. (2006). Each spectrum is an average of emission within the central 10" of each core. Each 
profile, asymmetrically blue with respect to the systemic velocity of each core (dotted lines), is indicative of inward motions, and can be 
reproduced well with a two-layer radiative transfer model (solid curves). 



Goodman (1998)of NH3 (1,1) emission over 3 isolated star- 
less cores, finding velocity gradients of ~0.3 to 1.40 km 
pc~^ with projected rotational axes not aligned along 
projected core axes. In addition. Swift et al. (2005) found 
a 1.2 km s~^ pc~^ gradient of NH3 emission aligned with 
the minor axis of the isolated starless core L1551-MC (see 
Figure 1). Furthermore, Lada et al. (2003) found evidence 
for differential rotation and angular momentum evolution in 
the B68 core, finding velocity gradients of 3.4 km s~^ pc~^ 
and 4.8 km s"! pc"! from C^^O (1-0) and N2H+ (1-0) 
emission respectively, which in turn likely trace outer and 
iimer core material respectively. 

Further investigations of velocity gradients include a 
study of 12 starless and 14 protostellar cores mapped in 
N2H+ (1-0) by Casein et al. (2002c), who derived a typical 
velocity gradient of 2 km s"^ pc^^ for both types of cores, 
and found that the spatial variations of the gradients indicate 
complex motions not consistent with simple solid body ro- 
tation. As found in the earlier study of cores by Goodman et 
al. (1993), these motions appear dynamically insignificant, 
since the derived rotational energies are at most a few per- 
cent of the respectively derived gravitational potential en- 
ergies. Even relatively small rotational energies, however, 
could influence protostellar formation during core collapse. 
For example, angular momentum evolution due to contrac- 
tion (see below) may induce disk formation or fragmenta- 
tion into multiple objects (e.g., see Hennebelle et al. 2004). 

In addition to plane-of-sky motions, molecular line pro- 
files can be used to trace the gas kinematics along the line- 
of-sight. For example, Walsh et al. (2004) determined that 
the relative motion between 42 isolated cores and their sur- 
rounding gas was relatively low (e.g., < 0.1 km s~^) by 
comparing profiles of a core tracer hke N2H+ with profiles 
of cloud tracers like C^^O and ^^CO. This result implies 
that cores do not move ballistically with respect to their sur- 
rounding gas. 



Molecular line profiles can also be used to detect inward 
motions in cores. Such motions can be detected using op- 
tically thick Unes, which will appear asymmetrically blue 
relative to symmetric, optically thin lines in a centrally con- 
tracting system. Indeed, inward motions may distinguish 
prestellar cores from starless ones, since they may indi- 
cate gravitational binding (see Section 2.1; Gregersen and 
Evans, 2000; Crapsi et al., 2005). Asymmetrically blue 
profiles were reproduced very well in collapse models of 
the protostellar envelope of B335 by Zhou et al. (1993), 
Choi et al. (1995), and Evans et al. (2005). In those pro- 
tostellar cases, the models indicated that the self absorption 
is due to the static envelope expected in the inside-out col- 
lapse {Shu, 1977) and the blue peak is stronger because of 
radiative transfer effects (see Zhou et al., 1993). Although 
such profiles may in principle also arise from rotational or 
outflow motions, 2D mapping and detailed modeling of the 
emission can distinguish between the different cases (e.g., 
Narayanan and Walker, 1998; Redman et al., 2004). In ad- 
dition, explanations other than inward motions would not 
favor blue-skewed over red-skewed profiles in a statistical 
sample. 

Figure 4 shows two examples of such profiles, seen in 
the Fi,F = 2,3 ^ 1,2 hyperfine component of the N2H+ 
(1-0) transition, observed towards the central positions of 
the L694-2 and LI 544 starless cores from Williams et al. 
(2006). Such a blue-skewed profile results from red-ward 
self-absorption, i.e., absorption along the line-of-sight by 
foreground, inward-moving gas of lower excitation. The 
statistics and the extended nature of many of the blue pro- 
files suggest that inward motions exist before formation of 
the central luminosity source, in contrast to predictions of 
the simplest inside-out collapse models (e.g., Shu, 1977). 

Such blue-skewed line profiles can be modeled with ra- 
diative transfer codes to estimate inward motion velocities. 
A simpUfied solution to the problem has been presented by 
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Myers et al. (1996) with a "two-layer" model, i.e., with two 
isothermal gas layers along the line-of-sight approaching 
each other, that can be used to determine rapidly the essen- 
tials of inward velocity. Figure 4 shows how well such mod- 
eling can reproduce observed asymmetrically blue profiles, 
with infall velocities of ~0.1 km s^^ for both cores. With 
the same goal, De Vries and Myers (2005) have recently 
provided a more realistic "hill" model, i.e., where the exci- 
tation increases with optical depth along the line-of-sight. 
Full radiative transfer solutions, on the other hand, can now 
be obtained relatively rapidly, and can be used to derive ve- 
locities (as well as densities, temperatures, and molecular 
abundances) of starless cores (e.g., see Keto et al., 2004; 
Lee etai, 2004b). 

Surveys of starless cores have revealed many candidates 
with profiles suggestive of inward motions. For example, 
Lee et al. (1999), using CS (2-1) as an optically thick tracer 
and N2H+ (1-0) as an optically thin counterpart in single- 
pointing observations, found 17 of 220 optically selected 
"starless" cores with asymmetrically blue profiles. Subse- 
quent modeling of these profiles suggested inward veloc- 
ities of 0.04-0.1 km s~^, i.e., subsonic. (The cores here 
are generally isolated and not found in turbulent molecular 
clouds.) Lee et al. (2004a) also observed 94 cores in the 
higher excitation CS (3-2) and DCO+ (2-1) Unes, finding 
18 strong inward motion candidates. The average inward 
velocity of 10 cores derived from CS (3-2) was found to 
be slightly larger than the average velocity derived from CS 
(2-1), i.e., 0.07 km s^^ vs. 0.04 km s^^ respectively, sug- 
gesting more dense, inner gas moves faster than less dense, 
outer gas in these cores, though the statistical sample was 
small. Indeed, CS, being depleted in the central and denser 
regions of the cores (see Section 4), only traces the external 
(radii > 7000 AU) envelope of the "chemically old" cores, 
so that the extended infall may reflect the (gravity driven) 
motion of the material surrounding these cores. 

The nature of the inward motions observed in starless 
cores can be probed in detail through multi-line analyses of 
individual objects. For example, the isolated core in L1544 
was found by Tafalla et al. (1998; see also Williams et al., 
1999) to have CS and H2CO line profiles suggesting inward 
motions over a 0.2 pc extent. This scale is inconsistent with 
that expected from thermal "inside-out" collapse of a sin- 
gular isothermal sphere since a central protostar of easily 
detectable luminosity would have hkely formed in the time 
required for such a collapse wave to have propagated so 
far (A CS/N2H+ mapping survey of 53 cores by Lee et 
al, 2001 also found 19 with infall that was very extended.) 
In addition, the inward velocities determined via two-layer 
modeling were up to ~0.1 km s^^. 

Although inward velocities derived from observations 
are larger than expected from some models of ambipolar 
diffusion, Ciolek and Basu (2000; CBOO) argued that such 
motions could indeed result from magnetically retarded 
contraction through a background magnetic field of rela- 
tively low strength, e.g., '^lO /iG, and tailored a specific 
ambipolar diffusion model for the L1544 core. Multiple 



lines of N2H+ and N2D+ observed from LI 544 by Caselli 
et al. (2002a) indeed showed some consistency with fea- 
tures of this model, including the velocity gradient observed 
across the minor axis and a slight broadening of line width 
and double-peaked line profiles at the core centre (possi- 
bly due to central depletion, but see Myers, 2005). In- 
deed, the ionization fraction of the LI 544 core found by 
Casein et al. (2002b) suggests its AD/FF timescale ratio 
« 1 (see Section 4.3). The observed fines at central lo- 
cations were wider than predicted from CBOO, however, 
suggesting possible modifications such as turbulence and, 
in particular, an initial quasi-static layer contraction due to 
Alfvenic turbulence dissipation (e.g., Myers and Zweibel, 
2001). Still more recent data have revealed further devia- 
tions from the CBOO model, however. For example, van der 
Tak et al. (2005) found evidence in double-peaked H2D+ 
1 10-1 11 line emission from the L1544 core for velocities 
(and densities) that are higher in its innermost regions than 
predicted by flie CBOO model. (The CBOO model could 
fit N2H+ and H^-''CO+ data of the core but they sampled 
larger radii than H2D+.) Such velocities could occur if the 
central region of the LI 544 core, where thermal pressure 
dominates, was smaller than predicted by CBOO model, or 
if collapse in these regions has become slightly more dy- 
namic (but still subsonic). 

Other isolated cores have shown inward velocities sim- 
ilar to those seen in L1544. For example, Tafalla et al. 
(2004) found evidence in N2H+ (1-0) Une emission toward 
the L 1 498 and L 1 5 1 7B cores for relatively high velocity gas 
(i.e., ~0.1 km s~^ offset from systemic velocities) that was 
coincident with regions of relatively high CS abundance. 
Such material may have been accreted onto the core only 
recently, suggesting a more stochastic growth of cores over 
time. In addition, Williams et al. (2006), comparing N2H+ 
1-0 of the L694-2 and L1544 cores, found similar inward 
velocity increases with smaller projected radii, although the 
L694-2 core has a shallower gradient but a similar central 
value to that of the LI 544 core (see Figure 4). 

Not all starless cores show kinematic behaviour similar 
to L1544, however. For example, Lada et al. (2003) found 
both blue- and red-skewed profiles of CS (2-1) across the 
B68 core grouped in patterns reminiscent of a low-order 
acoustic oscillation mode. Since the B68 core may be 
more in structural equilibrium than the L1544 core (e.g., 
see Alves et al., 2001), simultaneous inward and outward 
velocities of ~0.1 km s~^ magnitude in the B68 core may 
not have to do with core contraction and expansion. Instead, 
such oscillations may have been excited by external influ- 
ences, e.g., recent passage of a shock from expansion of the 
Loop I superbubble. 

Ambipolar diffusion models like that of CBOO may not 
be required to explain inward motions within starless cores. 
For example, Myers (2005) has recently proposed that the 
pressure-free, gravitational collapse of non-singular, cen- 
trally concentrated gas configurations alone is sufficient to 
explain observed velocity features of starless cores; i.e., no 
magnetic retardation is required. 
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5.2 lUrbulent Motions 

Pressure support against gravity within starless cores is 
provided by both thermal and turbulent motions of the gas, 
in addition to what magnetic field pressure can impart to 
neutral species through ions. Observationally, widths of 
fines from undepleted molecules in starless cores can be 
used to probe tfie relative influence of tfiermal and nonther- 
mal (turbulent) motions, and thus pressures in such cores, 
provided gas temperatures can be well determined (and 
magnetic field strengths; see chapter by Ward-Thompson et 
al.) Early NH3 observations showed that low-mass dense 
cores have typically subsonic levels of gas turbulence {My- 
ers, 1983). More recently, Goodman et al. (1998; see also 
Barranco and Goodman, 1998) found, in the cases they ex- 
amined, that the width of the NH3 lines within ^0. 1 pc of 
the core center does not follow the line width-size relations 
seen on large scales (e.g., oc BP'^; see Larson, 1981), 
suggesting that starless cores may constitute localized is- 
lands of "coherence." Goodman et al. further speculated 
that this reduction of turbulent motions may be the result 
of dense cores at 0. 1 pc providing a threshold for turbulent 
dissipation, e.g., through lower ion-neutral coupling from 
reduced internal ionization at high density and extinction. 

Regions of quiescent line width in other starless cores 
were analyzed by Caselli et al. (2002c), who found a lack 
of general correlation between the N2H+ (1-0) fine width 
and the core impact parameter in a sample of 26 cores, in 
agreement with earlier results. Furthermore, Tafalla et al. 
(2004; see also Tafalla et al., 2002) found a close to con- 
stant dependence between the NH3 or N2H+ line widths 
and impact parameter in the L1498 and L15 17B cores. This 
breakdown of the fine widtfi-size relation fias been repfi- 
cated in recent numerical MHD models of dense cores, 
described elsewhere in this volume; e.g., see chapter by 
Ballesteros-Paredes et al. Large numbers of the "cores" 
produced in these models, although even resembling BBS 
in density configuration, have larger internal velocities than 
the subsonic cores seen in Taurus. For example, Klessen et 
al. (2005) found only <23% of their model cores were sub- 
sonic. Comparisons between observed internal velocities of 
cores in other clouds and those predicted by turbulent mod- 
els would be very interesting. Regions of quiescence within 
clustered cores have also been noted by Williams and Myers 
(2000) and Di Francesco et al. (2004). (For further discus- 
sion of the evolution of turbulent motions in cores, see the 
following chapter by Ward-Thompson et al.) 

6. CONCLUSIONS 

This review described the numerous recent observations 
and resulting inferences about the nature of starless cores. 
These objects represent the state of molecular gas preced- 
ing gravitational collapse that defines the initial conditions 
of star formation. Recent observations, together with re- 
cent advances in theory and experimental chemistry, have 
shown that these objects reside at different points on evolu- 



tionary paths from ambient molecular material to centraUy 
concentrated, slowly contracting configurations with signif- 
icant molecular differentiation. Though sigruficant obser- 
vational strides have been made, sensitivity and resolution 
limits on past instrumentation have in turn limited our abil- 
ity to probe their internal physical and chemical properties. 
Fortunately, new instruments such as the wide-field fine and 
continuum focal-plane arrays on single-dish telescopes like 
the JCMT, LMT, and the TRAM 30 m, Nobeyama 45 m, 
and APEX Telescopes, and the new sensitive interferom- 
eters like the SMA, CARMA, and ALMA will be avail- 
able soon. These will provide the means to characterize 
numerous cores to high degrees, moving our understanding 
of their internal physical and chemical properties from an 
anecdotal to a statistical regime. 
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